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Abstract — Electronics components have an increasingly 
critical role in avionics systems and in the development of 
future aircraft systems. Prognostics of such components is 
becoming a very important research field as a result of the 
need to provide aircraft systems with system level health 
management information. This paper focuses on a 
prognostics application for electronics components within 
avionics systems, and in particular its application to an 
Isolated Gate Bipolar Transistor (IGBT). This application 
utilizes the remaining useful life prediction, accomplished 
by employing the particle filter framework, leveraging data 
from accelerated aging tests on IGBTs. These tests induced 
thermal-electrical overstresses by applying thermal cycling 
to the IGBT devices. In-situ state monitoring, including 
measurements of steady-state voltages and currents, 
electrical transients, and thermal transients are recorded and 
used as potential precursors of failure. 1,2 
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1. Introduction 

Predicting the Remaining Useful Life (RUL) of electronic 
components like IGBTs, MOSFETs (metal-oxide- 
semiconductor field-effect transistor), and BJTs (bipolar 
junction transistor) is one of the most challenging frontiers 
of Prognostic and Health Management systems. These 


components form the backbone of avionics systems that 
play an ever-increasing critical role in on-board, 
autonomous functions for vehicle controls, 
communications, navigation, and radar systems. 

During high-voltage operation, large internal electrical 
fields build up within power devices, which can increase the 
rate of degradation. For prognostic purposes the externally 
observable metrics of the device, such as collector 
current/voltage measurements during duty cycle on/off 
states of the IGBT that show degradation trends towards 
thermal runaway or latch-up, are correlated to physics-of- 
failure mechanisms like hot carrier injection, electro- 
migration, etc., and used in a particle filtering (PF) 
framework to carry out RUL prediction. 


2. Background 

A study of failure precursors for prognostics of IGBTs has 
been presented in [1]. This work focused on devices that 
had been aged under thermal overstress [3] resulting in 
identified precursors of failure that could be measured and 
computed from external signals. An accelerated aging, 
characterization and operational scenario simulation system 
was presented in [2]. This system allowed for aging under 
different operational conditions providing thermal and 
electrical overstress to devices, while monitoring in-situ key 
operational parameters that could serve as precursors of 
failure. The development of physics based models of aging 
and fault progression has also been considered. The work 
presented in [3] focuses on modeling the aging effects at a 
system level, by studying electrical power drivers using 
IGBTs. 

Failure Mechanisms Review 

Some major intrinsic faults relevant to transistor physics 
include dielectric breakdown, hot carrier injection, and 
electromigration [4, 5]. Some major extrinsic faults relevant 
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to transistor packaging-, include contact migration, wire 
lift, die solder degradation and package delamination [4, 6]. 

Dielectric breakdown occurs when a strong electric field 
induces a current channel through a previously insulated 
medium. Acute dielectric breakdown is typically the result 
of electrostatic discharge (ESD) and junction over-voltage. 
Time Dependent Dielectric Breakdown (TDDB) refers to 
the breakdown of gate oxide caused by chronic defect 
accumulation in the Si0 2 insulator. TDDB is shown to be 
advanced by increases in electric field strength [5]. A strong 
electric field may induce t energy into an electron or a hole 
which then becomes a “hot carried ’ due to the high kinetic 
energy stored [4]. Hot carriers have sufficient energy to 
tunnel and become trapped in gate oxide, and they are the 
primary cause of TDDB and contribute to device failures 
under normal operating conditions. Electromigration is a 
result of high current densities in silicon interconnects 
causing migration of metals. Formation of metal voids on 
interconnects can cause open circuits or high resistive paths, 
which in turn can result in poor performance or circuit 
malfunction [5]. 

Contact migration forms metal voids between external 
contact metals and the silicon. As metal voids grow, the 
aluminum or other metals can diffuse down to the silicon. 
This in turn can cause metal spikes to form deep in the 
silicon region which results in shorting the p-n junctions 
[7]. Wire lift occurs when the bond between the package 
wires connecting to the silicon die fail. Wire lift has been 
identified as a dominant failure mode in high power IGBTs 
[6]. Die solder degradation is another prominent package 
related fault. Solder attaching the silicon die and package 
heatsink develop cracks and voids due to thermal expansion 
mismatch between materials during expansion and 
contraction [6, 8]. 

Accelerated Aging Methodologies 

Thermal stress and electrical stress are the most common 
aging methodologies. Thermal cycling and chronic 
temperature overstress are prevalent thermal stress methods, 
with thermal cycling among the most prevalent accelerated 
aging methodology in electronics. Thermal cycling subjects 
devices to rapid changes in temperature differentials 
causing thermal expansion and contraction. Die solder 
degradation and wire lift are associated strongly with this 
aging method. Thermal overstress, another prevalent 
method, subjects devices to high temperatures for extended 
periods of time. TDDB is accelerated under high 
temperatures [9] and transistors have exhibited temperature 
dependant lifetimes accelerated by this mechanism [10]. 
IGBTs aged with self heating have shown changes in 
current ringing characteristics during switching [3]. 

Electrical overstress can be induced though transient and 
steady-state methods. Transient methods include electro 
static discharge (ESD), inductive switching and 


electromagnetic pulses. ESD is a leading cause of gate 
oxide failure [11] and hard switching of inductive loads, 
causeing voltage spikes which can cause significant damage 
to drain-source junctions [12]. Steady-state methods include 
chronic over-voltage and over-current. Applying high gate 
voltages [13], setting gate voltage ( V g ) to maximize drain 
current [14], and applying current overstress across the 
drain [15] have been shown to induce hot carrier and TBBD 

[9]. 


3. Prediction Algorithms 

Particle filtering is a sequential Monte Carlo (SMC) 
technique for implementing a recursive Bayesian filter using 
Monte Carlo simulations. It is primarily used for state 
estimation and tracking. The mathematical formulation for 
PF methods have been discussed in [16]. The basic idea is 
to develop a nonparametric representation of the system 
state probability density function (pdf) in the form of a set 
of particles with associated importance weights. The 
particles are sampled values from the unknown state space 
and the weights are the corresponding discrete probability 
masses. As the filter iterates the particles are propagated 
according to the system state transition model, while their 
weights are updated based upon the likelihood of the 
measurement given the particle values. Resampling of the 
particle distribution is done when needed in order to prevent 
the degeneracy of the weights. For state prediction purposes 
the same PF framework can be used by running only the 
model-based particle propagation step until the predicted 
state value crosses some predetermined end-of-life 
threshold. The predicted trajectory of each particle then 
generates an estimate of RUL, which can be combined with 
the associated weights to give the RUL pdf. The process is 
broken down into an offline learning part, and an online 
tracking and prediction part. During offline analysis, 
regression is performed to find representative ageing 
curves. Exponential growth models, as shown in Eqn. 1, are 
then fitted on these curves to identify the relevant decay 
parameters like C and 2: 

d=Cexp(-At) (1) 

where 0 is an internal model parameter of interest. More 
details of the PF framework used here can be found in [17]. 

The state and measurement equations that describe the 
semiconductor aging model are given below: 

Z k ^k - 1 • GXp { -Ajftk tfc-i ) } T CO k 

A* = A*_i + v k 

— \_Zk , A k\ 

y k = z k + v k (2) 

where the vector z, consists of the OFF state exponential 
time decay constants for the IGBT collector-emitter current 
(I C e), and matrices C and A contain their aging decay 
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parameters C and X values respectively. The z and A vectors 
are combined to form the state vector x. The measurement 
vector y comprises the time decay parameters inferred from 
measured data. The time index is denoted by k. The values 
of the C and A vectors learned from regression can be used 
to initialize the particle filter. The noise samples co , v and v 
are picked from zero mean Gaussian distributions whose 
standard deviations are derived from the given training data, 
thus accommodating for the sources of uncertainty in 
feature extraction, regression modeling and measurement. 
System importance resampling of the particles is carried out 
in each iteration, in order to reduce the degeneracy of 
particle weights. This helps in maintaining track of the state 
vector, even under the presence of disruptive effects like 
unmodeled operational conditions (in our case, high 
temperatures). 

The system description model developed in the offline 
process is fed into the online process where the particle 
filtering prognosis framework is triggered by a diagnostic 
routine. The algorithm incorporates the model parameter as 
an additional component of the state vector and thus, 



Figure 2 - IGBT Voltage and Current Waveforms 

At high temperature thermal cycling these current decay 
curves change due to hot carrier injection. The trends in 
these curve shifts can act as a precursor to catastrophic 
latch-up and subsequent thermal runaway conditions. Figure 
3(a)-(c) shows the shift in the OFF-state current decay as 
aging progresses, finally ending in a latch-up condition as 
shown in Figure 3(d). 
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4. Experiment Details 


The IGBT behaves essentially like a MOSFET and a BJT 
connected in a Darlington configuration. The input of the 
IGBT is modeled by the MOSFET, while the output is 
modeled by a BJT. As a combination of these two devices, 
the IGBT takes on their characteristics as well as some 
additional ones caused by their interaction. For example, it 
exhibits a tail collector-emitter current (I C e) while turning 
off, which is modeled by a decaying current flowing 
through the BJT after the MOSFET has been turned off. 
This is shown in Figure 2. 
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5. Results 
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Figure 3 - (a), (b) and (c) Shift in IGBT Off-state 
Current Decay with Aging; (d) Latch-up Condition 

Accelerated Aging Settings 

Thermal overstress accelerated aging was conducted on 
IGBTs using the agile accelerated aging, characterization 
and scenario simulation system for gate controlled power 
transistors presented in [2]. An International Rectifier 
IRG4BC30KD IGBT with a 600V/15A rating in a TO220 
package was aged with no external heatsink. The collector- 
emitter terminals were connected in series with a load 
power supply and a 0.2Q load resistor. An infrared sensor 
was used for temperature measurement of the IGBT case. 
The gate signal was chosen to be 10V DC. The aging 
process was controlled by a hysteresis temperature 
controller with set points of 268°C (low) and 270°C (high), 
switching the gate voltage ON and OFF as a control action. 
The load power supply voltage was 1 0 V which resulted in a 
current of 8 A. An additional temperature threshold 
controller, with a set point of 305°C, was programmed to 
turn OFF the load power supply and end the experiment in 
the event of thermal runaway and latching failures. The 
IGBT under consideration was aged for -210 minutes until 
latch-up resulting on thermal runaway. 


During the model learning process, we extract the tail 
current sections from the I C e waveform, transform them to 
the log domain and then fit 3 rd degree polynomials, as 
expressed below: 

Ice (0 “ expfc/ + P 2 t * 2 * * + P 3 t + P 4 }. (3) 


The equivalent fits in the time domain are denoted by the 
cyan lines in Figure 3(a)-(c). The fit parameters P b P 2 , P 3 
and P 4 are then tracked over the entire aging duration. It is 
to be noted that the parameters are highly correlated with 
each other (as denoted by the correlation coefficient matrix 
shown below in Eqn. 4), and consequently, we use only Pi 
as the parameter to track. 


1.0000 - 0.9964 0.9873 - 0.9747 

- 0.9964 1.0000 - 0.9972 0.9899 

0.9873 - 0.9972 1.0000 - 0.9978 

- 0.9747 0.9899 - 0.9978 1.0000 


(4) 
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Figure 4 - Trending of Curve Fit Parameter Px 

Overall, the parameter shows an exponential growth rate 
(negative) as shown in Figure 4, possibly indicating an 
Arrhenius aging process. Near the end, the curve shows 
some anomalous behavior as the IGBT approaches the 
latch-up condition. Originally the Arrhenius equation 
described the temperature dependence of the rate of a 
chemical reaction. Currently it is perceived largely as a 
empirical relation which can be used to model the 
temperature variance of many thermally-induced processes. 

In running the particle filter framework, we used an 
exponential growth model for Pi like the one shown in Eqn. 

2. Since we did not have separate learning and testing 

datasets, we did not use the regressed C and X values from 

the model learning process to jump start our PF algorithm. 

Instead, we use 0 and 1 as the initial values for Pi and / P| . 
The EOL threshold is chosen as -2.5x1 0 5 , which is 
approximately the value of Pi at the end of its exponential 
behavior. Without loss of generality, the prediction time is 
arbitrarily chosen to be 51.875 minutes. 
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Figure 5 shows the performance of the PF in tracking and 
prediction modes. The inset plot shows the 2 Pi parameter 
estimation. After the prediction timeline, the mean of the 
last 100 2, P values is used for the propagation of each 
particle according to the state transition model (first line in 
equation 3). No additional computed Pi values are used to 
update the particle weights. The propagated particle values 
are compared against the EOL threshold to compute the 
distribution of time at EOL. A mixture of Gaussians is then 
fitted in a least square sense to these EOL values to plot the 
magenta EOL pdf. The RUL pdf is simply the prediction 
time subtracted from the EOL pdf. It is to be noted that 
despite the generic initial particle values, modeling 
uncertainties and system noise, the PF performance is very 
good. 

6. Conclusions 

This prognostic exercise, although very preliminary, gives 
rise to the hope that one can reliably detect failure 
precursors in semiconductor device behavior and use them 
in an intelligent prediction framework (e.g. particle 
filtering) to derive RUL estimates. The fact that the PF 
provides a RUL pdf, instead of a single MTBF value, makes 
the interpretation of the prognostic result more useful and 
intuitive. 

There are some significant limiting assumptions employed 
here: Simple exponential growth models wree used to 
explain the aging behavior of IGBTs. Although the 
prediction results are very good, such an approach has 
limitations when applying the aging parameters and EOL 


thresholds learnt from one device to another. In order to 
achieve that, one needs to explicitly derive physics-of- 
failure based aging models, identify their parameters from 
externally observable device characteristics and then use 
those models in a PF framework to carry out prognostics. In 
such a case, aging behavior will be linked to the changes in 
the internal model parameters which can then be tied to 
failure mechanisms like hot carrier injection, electro- 
migration etc. The focus of our future work will be to 
devise aging experiments to fulfill the above objectives. 
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